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An electrical circuit consisting of two capacitively coupled inductive loops, each interrupted by a 
Josephson junction, is analyzed through the classical RCSJ model. The same circuit has recently 
been studied experimentally and the results were used to demonstrate quantum mechanical en- 
tanglement in the system by observing the correlated states of the two inductive loops after initial 
microwave perturbations. Our classical analysis shows that the observed phenomenon exists entirely 
within the classical RCSJ model, and we provide a detailed intuitive description of the transient 
dynamics responsible for the observations. 



I. INTRODUCTION 

In 2005, McDermott et al. [l[ reported the experimen- 
tal observation of quantum mechanical entanglement in a 
system of two qubits, each comprised of a superconduct- 
ing loop interrupted by a single underdamped Josephson 
junction, and weakly coupled via a capacitor. A sec- 
ond report with some additional detail soon followed Q, 
and a review by Siddiqi and Clarke of the key ideas 
appeared in that same issue of Science. 

Briefly, the model put forward to explain the exper- 
imental observations was as follows. The potential en- 
ergy of a superconducting ring interrupted by a Joseph- 
son junction, plotted as a function of the junction phase 
variable, can possess a shallow well if biased with an ap- 
propriate magnetic field. Quantum mechanics dictates 
discrete energy levels within a finite well, and a suf- 
ficiently shallow well would have just a few permitted 
bound states, of which the two lowest are denoted |0), 
The argument has been advanced that such a sys- 
tem can be viewed as a fictitious particle being limited 
to these quantum states [H, rather than, as in a "classi- 
cal" picture, freely exploring the potential surface. This 
distinction is crucial to what follows. 

In the experiments, two loops were subjected to a 
steady flux bias that made the occupied wells shallow. 
One of the loops was given additional time-dependent 
flux biasing in the form of a 10 ns microwave burst. The 
microwave frequency was set to match the separation of 
the two levels within the driven well, so it was imag- 
ined that the pulse would selectively populate the first 
excited state of the driven loop, initially putting the cou- 
pled system in the 1 01} state. Then, after a "free evo- 
lution time" during which the qubits interact, the states 
of both loops were probed using identical measurement 
pulses of exactly the correct amplitude to reduce the bar- 
rier and facilitate selective tunneling (or, more precisely, 
Macroscopic Quantum Tunneling (MQT)) out of the up- 
per state |1), but not out of the lower state |0). Hence, if 
either well were in an excited state, then escape via tun- 
neling would lead to time dependent phase dynamics and 
a resulting observed signal; if it were in the lower state, 
no signal would be seen. This probing was repeated many 



times for each selected tf ree and the accumulated data 
were argued to reveal the probabilities for |00), 1 01) , |10), 
and 1 11), and to confirm that entanglement had been ob- 
served. 

We here provide a classical Resistively and Capaci- 
tively Shunted Junction (RCSJ) analysis of the capac- 
itively coupled qubit system. This builds upon our pre- 
vious work which has reproduced several experimental 
observations, such as those seen in y|, including: multi- 
peaked distributions [||, Rabi-oscillations Q, Ramsey 
fringes [1] , and spin-echo @ . In this model the microwave 
pulse stimulates small phase oscillations in the driven 
loop and, because of the coupling, the passive loop soon 
develops its own phase oscillations. Thus we study a 
straightforward picture - that of a pair of weakly cou- 
pled oscillators. Our analysis shows that these coupled 
phase oscillations can yield very good quantitative agree- 
ment with the core of experimental observations. It also 
points to the crucial role of thermal noise in this system. 



II. THE MODEL 

We begin with the circuit shown in Fig. 1. Each junc- 
tion is characterized by a critical current I c , resistance R, 
and capacitance C; each loop has an inductance L. The 
two loops are coupled through a capacitance Cx- Let 
&i x be the externally applied flux on Loop #1 and $>2x be 
the externally applied flux on Loop #2. With overdots 
denoting derivatives in dimensionless time r = wji, and 
with the junction plasma frequency wj = y f 2eI c /HC ', the 
resulting equations of motion for the two junction phases 
can be shown to be 

ipi + aifii + sin<£i = J x (tfi2 - fti) - Pl 1 (<fi + 2%Mi x ) 

(1) 

ip2 + o>ip 2 + sin ip 2 = jxifii - 02 ) - Pi 1 (<P2 + 2nM 2x ) 

(2) 

where a = 1/uijCR, and 0l = 27rL/ c /$o; $o being the 
flux quantum. Mi w — <&i X /&o is the normalized applied 
flux through loop i and the mutual coupling coefficient 
is J x — C x /C. The characteristic energy of this system 
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FIG. 1: Configuration of two capacitively coupled supercon- 
ducting loops, each containing a Josephson junction repre- 
sented by the RCSJ equivalent circuit and an associated noise 
source. 



is Ej = J c $o/27T. The equations of motion can also be 
put in the convenient form, 

tp a + aip a + sin ip a cos ipb = -fil 1 [ip a + 2irM a ] (3) 
g~ x (pb + oi^Pb + sin ip b cos (p a = 1 [fb + 2nM b ] , (4) 

where ip a = (ipi + <P2)/2 and ipb = (ipi — f>2)/2 are 
transformed variables, and M a = {M\ x + M2 X )/2, Mb = 
(Mi x — M2x)/2, an d g^ 1 = 1 + 2~f x are the corresponding- 
magnetic fields and coupling, respectively. Based on pub- 
lished data, we set the parameters at a — 5 x 10 -5 (very 
light damping), f3 L = 2.841, and g = 0.9954, I c = 1.1 //A, 
C = 1.3 pF, uj^ 1 = 0.02 ns. Both loops are biased with a 
dc flux 0.6941 (resonance «5.1GHz) and with superim- 
posed pulses as shown in the upper panel of Fig. 2. 

III. SIMULATION RESULTS 

We have numerically solved the coupled equations ([3]) 
and ([4]) using both fourth order Runge-Kutta and Verlet 
algorithms. Simulations at zero temperature were con- 
ducted as follows. 

The duration of the microwave burst was set at 100 
plasma periods; the amplitude and normalized frequency 
were 0.000580 and 0.989, respectively. Our results are 
shown in the three lower panels of Fig. 2. In each plot, 
the darker trace is the phase in Loop #2. The effect 
of the microwave burst is easily seen - the phase <f2 is 
kicked into rapid oscillations around its initial rest value 
of — 1.52 rad. These oscillations continue after the mi- 
crowave burst has ended. Because of the coupling, the 
phase in Loop #1 also oscillates. Furthermore, in a 
manner characteristic of weakly coupled oscillators [Io| , 
each waveform has the appearance of an amplitude mod- 
ulated carrier, and the envelope of one is exactly out of 
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FIG. 2: (a) Bias sequence for the two loops (not to scale). 
(b)-(d) Numerical solutions of the coupled loop equations. 
In each plot, the darker wave form is for Loop #2, the lighter 
for Loop #1. In (c) and (d), a rotation pulse is applied at 
time 9550. The insets show portions of the upper envelopes 
of the modulated phase oscillations. 



phase with the other. Note that the amplitude of these 
modulations is very small compared to the dc level. The 
junction phase oscillations are so rapid (the frequency ra- 
tio of carrier to modulation is about 400: 1) that on the 
scale of this figure they become compressed into solid 
shading. 

A rotation pulse of amplitude 0.000915 was applied at 
time 9550. For these simulations, we chose widths of 0, 
60, and 120 plasma periods, equivalent to ns, 7.5 ns, and 
15 ns. The insets in each of the three panels show the up- 
per modulation envelopes as functions of time and their 
responses to the different rotation pulse widths. The sig- 
nificance of these plots is their exact matching to the time 
dependence of the experimental probabilities Poi an d Pio 
given in Fig. 2B,C,D of Q. Steffen et al. say: "The oc- 
cupation probabilities Poi and Pw oscillate out of phase 
with a period of 100 ns, consistent with the spectroscopic 
measurements." Our numerical results exhibit the same 
shapes and phase relationships for the envelopes. As well, 
our simulation data show the period of the envelopes to 
be 87ns. In other words, the modulation envelopes that 
are a natural outcome of the classical physics of weakly 
coupled oscillators are a match to the experimental data, 
and the underlying dynamical phenomena of both the 
modulation and the results of manipulation observed in 
Ref. Q are well represented in the RCSJ model. 

The experiments reported in were conducted at 25 
mK . Following their experimental protocols, we added 
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FIG. 3: Numerical simulations for T = 25 mif with thermal- 



ization time of 10 time units, and a =5x10 



Pulsed mi- 



crowave frequency is 0.625, and pulse duration and amplitude 
are 507 and 0.00058, respectively. Each marker represents an 
average of 1,000 measurements conducted with a triangular 
probe pulse of amplitude 0.0538 and width 500. (b) and (c) 
rotation pulse application time is 818.75 after termination of 
the microwave, and rotation amplitude is 0.000915. (b) 90° 
rotation pulse duration is 253.5. (c) 180° rotation pulse du- 
ration is 507. 
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FIG. 4: Analyses of system modes with parameters given in 
the text. Single loop resonance (a) and energy (b) as a func- 
tion of oscillation amplitude A. (c) Modulation frequencies; 
simulated (markers) and perturbation analysis (lines). Dif- 
ferent modes are engineered through the initial conditions 
i^i(O) = ipo + Ai such that mode-I is obtained by A\ — A 
(small) and A2 = 0, mode-II by A\ = A (large) and A2 = 0, 
and mode-Ill by Ai ~ A2 (large) . Examples of the three char- 
acteristic modes observed in (c): mode-I (d), mode-II (e), and 
mode-Ill (f). 



will discuss after having identified the dynamical modes 
in the system. 



appropriate noise to the system equations © and Q, 
simulating 25mK thermal noise with a short transient 
time of 1,000 time units. For any chosen moment fol- 
lowing the microwave burst (tf ree ), repeated simulation 
runs were carried out and data were gathered on how 
often the individual loops experienced an escape follow- 
ing the application of the probe pulse (amplitude 0.0358). 
The results are shown in Fig. 3. The close correspondence 
with the reported experimental results in [2| is obvious. 

Comparing Fig. 2 and Fig. 3, it is apparent that in com- 
bination with noise, the probe pulse effectively 'teases 
out' the modulation profiles of the junction phase oscil- 
lations. Equivalently, the likelihood that the "particle" 
will escape from its well will be higher at moments when 
the noise-free oscillation amplitude (energy) is relatively 
large, and conversely, when the oscillation amplitude is 
small, the chance of escape is less. Hence the classical de- 
scription of the system, embodied in Eqs. (|3l4p captures 
all the significant features observed in the experiments. 
However, there is a wrinkle in the above picture that we 



IV. PERTURBATION ANALYSIS 

Let us first provide a simple analysis for the single 
qubit. Inserting the ansatz (p = ipo + A sin cut into Eq. ([1]) 
for "fx — we get the following relationships between the 
fixed phase ipo and the resonance frequency u r to the 
magnetic field M x and amplitude A of oscillation: 

^ = PI 1 + {J (A) + J 2 {A)) COS (5) 

where 

J (A)f3 L sin^o = -<Po- 2ttM x . 

The latter expression provides the mean phase (po, while 
the former Eq.© is the anharmonic resonance shown in 
Fig. 4a for small to moderate oscillation amplitudes. 

With this parameterized ansatz, we can provide the 
energy Hr, of a single qubit as the averaged sum of kinetic 
energy and potential energy U(ip) = 1 — cos tp + [tp + 



4 



2irM x } 2 /2(] L 



H L = (-^ + U^)) = -AW+ 

-^—A 2 + 1 - Jo (A) cos <A) + Jh^o + 2nM x } 2 



(6) 



This energy is shown in Fig. 4b for the upper of the two 
potential wells of the qubit potential. Notice that the 
energy of the fixed-point is H = U(ipo). With these 
basic anharmonic relationships we can now proceed to 
analyzing the inherent system resonances. 

Mode-I: beat frequency mode. In Eqs. ([3]) and |4}, 
with Mb = 0, a ~ 0, and small amplitude oscillations of 
ipi and if2, we write (p a — (po+ipa, where ipo is a constant 
and \ip a \ 1, and we will similarly assume |c/?f,| <C 1. 
Inserting this ansatz with a = (for simplicity) gives 
the linear equations 

i) a + sin Lp + (cos + (3Z l )ij)a + PI 1 (<fo + 27rAf a ) = 

(7) 

g~ x (pb + (cos ip + Pl 1 )^ = . 

(8) 

From this we directly obtain [3l sin^o + <^o + 2nM a = 0, 
and we get the linear resonance frequencies uj a and u>b 



Mode-III: phase locked mode. This is a mode, il- 
lustrated in Fig.4f, which can be analyzed similarly to 
the classical Rabi-type oscillation outlined in Ref. 0]. 
We define the ansatz tpi = ipo(Ai) + Ai sm(ui r (Ai)t + 0i) 
(i = 1, 2), A\ Ki A2 w A, and 9\ — —62- Considering the 
energy flow into qubit #1 from the coupling, 

Hi = 7^1(^2 - <f>i) (11) 
ps Aui r cos(ui r t + 81) (12) 
x [Auj 2 sin(uj r t + 61) - Aco 2 sin(uvi -61)] , 

we write the energy Hi of qubit =f/=l using Eq. The 
energy change AHi of qubit #1 over one time unit can 
therefore be expressed 

AiJr = ^ = d ^^0\ 



dt 
1 



dAi duj r 



7/, = ^Af^sm^) (13) 



dH x dA •• 2 3o 



(14) 



dAi dto r 

for \9i I <C 1. This provides a slow modulation frequency 



(15) 



of ^ a and ^, respectively, as w Q = cos^o + L and where we nQte that duj/dA < Q 



gtu 2 . Thus, in addition to the microwave carrier 



frequency ss uj a , we see a modulation frequency at 



■ Wft = (1 - y^O 



(9) 



which, with the reported experimental parameters, is 
given by uj a « 0.63199 and 0/ w 0.001455. Notice 
that 2tv/Hj is the exact modulation period (86.4ns) 
observed in both experiments and simulations. The 
energy modulation of this mode is exemplified in Fig.4d. 

Mode-II: phase-slip mode. A direct consequence of 
the anharmonicity of the system is that mode-I can only 
exist for low system energy. The reason is that the 
two different oscillation amplitudes {A\,A2) observed in 
Fig.4d give rise to two different resonance frequencies 
seen in Eq. ([5]) and Fig. 4a. Thus, the linear mode-I 
cannot exist if the disparity between the resonance fre- 
quencies overwhelms the linear modulation frequency Qi. 
If that happens, then the two oscillators will shift to a 
mode-II, where the two loops will advance in two differ- 
ent energy and frequency states, as illustrated in Fig.4e. 
The observed energy modulation is a result of the mutual 
phase-slip between the oscillators, which for systems of 
different oscillation amplitudes A\ > A2 can be approxi- 
mated by 



n 
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u r (A 2 ) - Lo r {Ai) 



(10) 



The above perturbation analysis provides definite ex- 
pectations for the type of excitations one can expect and 
how they relate to the experimental observations. Judg- 
ing from Figs.4d-f, we infer that only mode-I can be the 
experimentally observed, due to the statements in Ref. Q 
about using identical probe pulses to get the presented 
high resolution and fidelity. The complexity of the mod- 
ulation resonances are shown in Fig. 4c with comparisons 
between evaluated modulation frequencies and direct nu- 
merical simulations for a — as a function of oscillation 
amplitudes. We see that mode-I is always present for 
small A, and that it is never present for large A due the 
anharmonic modes II and III. The transition to the an- 
harmonic modes appears at a system energy given by the 
the magnitude of A - see Figs. 4c and 4b, where we have 
also indicated the energy level of 25mK. It is here appar- 
ent that mode-I does not exist with energy content much 
larger than the thermal energy for the given system pa- 
rameters. We have numerically verified this for a fully 
thermalized system at 25mK. Only if not fully thermal- 
ized (as is the case for Fig. 3 above) or at thermodynamic 
temperatures below 3mK do we observe mode-I. 



V. DISCUSSION 

We note that results of the classical model depend on 
the choice of the phenomenological damping parameter 
a, which is not directly given by the experimental data 
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in [2|. The inverse of this parameter relates to both the 
decay time for coherent signals in the system and the 
characteristic time of thermalization; the value of this 
parameter is therefore a significant component to un- 
derstanding the system behavior. For this presentation 
we have chosen a very small value of damping in line 
with the observed decay times of the modulations as well 
as with the assumed large sub-gap resistance associated 
with high quality aluminum junctions at the very low 
temperature used for the experiments. 

With the long history of the RCSJ model explaining 
experimental observations, as outlined in the introduc- 
tion, it is natural to approach any new configuration from 
the same starting point. For the present system, such di- 
rect classical analysis provides excellent agreement with 



what is actually seen in the experiments, as illustrated in 
Figs. 2 and 3. However, our analysis also points to a still 
unresolved question of thermal effects and the stability of 
the relevant modes in the balance between temperature 
and anharmonicity. 
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